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Abstract. The gonadal cycles of four groups of eastern 
oysters, Crassostrea virginica (Gmelin), including native 
stocks collected that year and inbred strains (reared in 
Delaware Bay for 5-6 generations) from both Long Island 
Sound and Delaware Bay, were compared in Delaware 
Bay in 1987. Inbred strains resembled their respective na¬ 
tive stocks; both Long Island groups initiated gonadal de¬ 
velopment and spawning about one month earlier and 
spawned over a shorter duration than both Delaware Bay 
groups. Analysis of covariance revealed that the effect of 
time on gonadal development was statistically different 
(P < 0.05) for all between-location group comparisons, 
but not for the two within-location comparisons. Thus, 
after six generations of inbreeding in Delaware Bay, Long 
Island oysters maintained their characteristic pattern of 
gonadal development and spawning, indicating the exis¬ 
tence of genetically different environmental requirements 
for gonadal maturation between the two locations. 

Introduction 

In t rasped tic variation in various aspects of reproduc¬ 
tion has been noted for several species of marine bivalves, 
and can be either genetic or non-genetic (adaptive) in na¬ 
ture (see Sastry, 1979; Barber and Blake, 1991). The sug¬ 
gestion that there are genetically distinct populations or 
“races"' of eastern oysters, Crassostrea virginica (Gmelin), 
was first made by Loosanoff and Engle (1942). Stauber 
(1950) concluded that there were three physiological races 
of oysters, based on differences in temperatures at which 
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spawning was initiated in Long Island, New Jersey, and 
Texas populations. Because transplant experiments were 
not conducted, the relative contribution of genetic versus 
non-genetic factors to the observed differences was not 
determined. Loosanoff and Nomejko (1951) and Loos¬ 
anoff (1969) transplanted stocks from various locations 
along the eastern and Gulf coasts of the U.S. to Milford, 
Connecticut, and histologically compared gonadal devel¬ 
opment. Differences in the timing and extent of gonadal 
development and spawning between some of these stocks 
in the common environment suggested the existence of ge¬ 
netic differences in the environmental factors regulating 
gonadal maturation. The possibility of acclimatization 
after transfer to the new environment, however, was not 
precluded. As stated by Loosanoff (1969), the best way to 
separate the genetic and environmental aspects of oyster 
reproduction would be to use “successive generations of 
laboratory-reared oysters originating from parents of dif¬ 
ferent geographical areas." 

The oyster (C. virginica) selective breeding program of 
the Rutgers University Shellfish Research Laboratory in¬ 
cludes a variety of wild (imported) stocks and laboratory- 
reared strains that originated from several locations (Has- 
kin and Ford, 1979; Ford and Haskin, 1987). Initial ex¬ 
amination of some of these strains (originating from Long 
Island Sound, Delaware Bay, and James River, Virginia, 
stocks) indicated that after several generations of inbreed¬ 
ing and maintenance in Delaware Bay, the timing of spe¬ 
cific reproductive events remained characteristic of the 
site of geographic origin (Ford et ai, 1990). The object 
of the present study was to establish whether these ob¬ 
served differences are attributible to intraspecitic genetic 
variation. This was accomplished by simultaneously 
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comparing cycles ot gonadal maturation and spawning in 
both wild stocks and inbred strains originating from two 
locations: Long Island Sound and Delaware Bay. 

Materials and Methods 

Native oysters from Long Island Sound (LIN), native 
oysters from Delaware Bay (DBN), a 6th generation, 1985 
year class, inbred strain originally imported from Long 
Island Sound to Delaware Bay in 1964 (Lll), and a 5th 
generation, 1985 year class, inbred strain originating from 
Delaware Bay (DBI) were obtained in March 1987. Be¬ 
cause the endoparasite Haplosporidiwn nelsoni (MSX) 
alters gonad development and relative fecundity (Barber 
et al, 1988; Ford and Figueras, 1988), all oysters used in 
this study were first “purged” of MSX (Ford and Haskin, 
1988a, b) by being suspended in wire trays from the lab¬ 
oratory dock (Bivalve, New Jersey) in the Maurice River 
(salinity < 5 ppt) for several weeks. On 1 May 1987, all 
experimental groups were moved to the tidal flats of lower 
Delaware Bay, Cape May County (salinity 18-24 ppt), 
for the remainder of the study. 

Each group of oysters was sampled periodically between 
17 March and 5 October 1987; ten individuals were 
shucked and fixed in Davidson's solution. A standard 
transverse (anterior) section across gill, stomach, intestine, 
and digestive diverticula was dehydrated, cleared, and 
embedded in paraplast. Six-pm sections were stained with 
iron hematoxylin, acid fuchsin, and aniline blue. Slides 
were examined for MSX prevalence (% of sample infected) 
and intensity (Ford, 1985; Barber et al, 1988). Gonadal 
development was assessed from the histological sections 
using the Bioquant Image Analysis System (R&M Bio¬ 
metrics, Inc., Nashville, Tennessee) to obtain a Gonadal 
Area Index, defined as the ratio of gonadal area to the 
entire visceral mass area X 100. The gonadal area tech¬ 
nique is a sensitive indicator of gametogenic events in 
oysters (Mori, 1979; Barber et al, 1988) and is more pre¬ 
cise than measuring the thickness of the gonadal layer 
alone, as performed by Loosanoff and Nomejko (1951). 

Cycles of gonadal development of the four groups of 
oysters, based on mean gonadal area indices, were com¬ 
pared using analysis of covariance after correcting for het- 
eroscedasticity with the least squares method ol White 
(1980). Analysis of covariance was accomplished using 
the dummy variable regression procedure of Wonnacott 
and Wonnacott (1970). Gonadal area index was regressed 
on a quadratic function of time (elapsed days from the 
date of the first sample) as: 

Gonadal Area Index 

= c v. Days 4- ofjDj Days + ft Days~ T ftD 3 Days**, 

using the Tobit estimation to overcome the problem of 
zeros in the data set (Goldberger, 1964). The t statistic 


was then used to test the null hypothesis that and 

ft - ft, for each pair-wise group comparison. Conven¬ 
tional analysis of variance was not applied because of the 
presence of zero-valued observations, which invalidated 
the required normality assumption. Moreover, conven¬ 
tional analysis of variance compares only means, while 
the procedure used here allowed a global comparison. 

Water temperature data from the experimental site were 
collected only sporadically during 1987. Weekly averages 
of “calculated” daily water temperatures were therefore 
obtained from the regression of daily air temperatures 
(Cape May, New Jersey) on actual water temperatures 
collected at the experimental site from 1978-80 and 1985- 
88 (n - 358; r = 0.77). These are presented in Figure 1. 

Results 

The haplosporidan parasite II. nelsoni had little effect 
on gonadal development in any of the groups in this study. 
Prevalence and intensity of the parasite remained rela¬ 
tively low throughout the first four months of the study 
period because of the exposure of oysters to low salinity 
prior to initial sampling (Table 1). Prevalence of H. nelsoni 
did not exceed 50% until 15 July, and systemic infections 
exceeded 50% only in the LIN group, after 15 July and 
the completion of spawning in that group. Gonadal de¬ 
velopment in C. virginica is most affected in individuals 
having systemic II. nelsoni infections (Barber et al, 1988; 
Ford and Figueras, 1988; Ford et al., 1990). 

Cycles of gonadal maturation and spawning for the four 
groups of oysters are represented as mean gonadal area 
indices in Figure 2. Both Long Island groups (LIN and 
LII) initiated gonadal development in April and had 
maximal gonadal areas in late May (Figs. 2A, B). Some 
spawning occurred in both Long Island groups in late 
May and June, but greatest spawning activity occurred in 
July. Long Island oysters (both LIN and LII) were spent 
by August. Although some gonadal development was ap¬ 
parent in the DBI group in late April, most gonadal growth 
in both Delaware Bay groups (DBN and DBI) occurred 
in May (Figs. 2C, D). Maximal gonadal areas were found 
for both Delaware Bay groups in June, which was followed 
by a period (from June to August) of what appeared to 
be partial spawning and redevelopment. Final spawning 
and the cessation of gonadal activity in the Delaware Bay 
groups occurred in September. Thus there were similarities 
between inbred strains and native stocks within each site 
of origin but differences between sites of origin. 

The null hypothesis of similarity of coefficients was re¬ 
jected in all four pair-wise comparisons between strains 
having different geographic origins (Table II), indicating 
that the timing of gonadal growth and spawning of both 
native stocks and inbred strains from Delaware Bay were 
statistically different (P < 0.05) from those of oysters (both 
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1987 

Figure 1 . Weekly mean vvaler temperature at lhe experimental site, lower Delaware Bay, as calculated 
from the regression of daily air temperature at Cape May, New Jersey, on water temperature readings from 
1978 to 1980 and 1985 to 1988. 


native and inbred) from Long Island Sound. On the other 
hand, the null hypothesis was accepted in the two com¬ 
parisons involving inbred and native groups having the 
same site of origin (LIN-LI1 and DBN-DB1), indicating 
that the timing of gonadal maturation cycles was similar 
for native stocks and inbred strains within each site of 
origin. 

Discussion 

There were distinct differences in the timing of gonadal 
development and spawning of both native stocks and 
inbred strains between the Long Island and Delaware Bay 
sites of origin. Oysters in the Long Island groups initiated 
gonadal growth, achieved maximal gonadal development, 
and began spawning about one month earlier than oysters 
in the Delaware Bay groups. Additionally, both Delaware 
Bay groups exhibited a protracted period (about three 
months) of partial spawning and redevelopment. The sta¬ 
tistical comparison of gonadal area indices reinforced the 
observations reported in this study and those of Ford el 
al. (1990). Time (elapsed days) had the same effect on 
gonadal area index for the LIN and LI1 groups and for 


the DBN and DBI groups, but not for any between-lo- 
cation comparison. Thus, the timing of gonadal devel¬ 
opment differed significantly between the Long Island and 
Delaware Bay locations. These differences were main¬ 
tained in the Long Island inbred strain, even after six 
generations (23 years) in Delaware Bay, demonstrating 
that there are genetically determined differences in en¬ 
vironmental criteria necessary for the initiation or com¬ 
pletion of a particular gametogenic event at these two 
locations. 

Several potential environmental regulators of bivalve 
reproduction have been identified, with the foremost of 
these being temperature (Sastry, 1979; Barber and Blake, 
1991). Scallops (Argopeclen irradiatis) from Woods Hole, 
Massachusetts, and Beaufort, North Carolina, acclimated 
to similar laboratory conditions prior to the initiation of 
gonadal development, exhibited maximal gonadal growth 
at 15 and 23°C, respectively, indicating a non-adaptive 
(genetic) difference in temperature requirements for go¬ 
nadal development between these two populations (Sastry, 
1966). One of the original criteria upon which the deter¬ 
mination of physiological races of oysters was based was 
the temperature at which spawning began. According to 
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Table I 


Prevalence (P, %) and intensity (I, expressed as % systemic infections) 
o/Haplosporidium nclsoni in experimental groups of oysters, 
Crassostrea virginica, on the various sampling dates 


Date (1987) 

LIN 

%P/%J 

LI1 

%P/%I 

DBN 

%P/%I 

DBI 

%P/%1 

17 March 

_ 

20/0 

50/20 

0/0 

29 April 

0/0 

10/10 

20/0 

0/0 

7 May 

0/0 

0/0 

— 

— 

] 3 May 

0/0 

0/0 

0/0 

0/0 

21 May 

0/0 

0/0 

0/0 

0/0 

27 May 

0/0 

20/10 

0/0 

0/0 

3 June 

0/0 

20/10 

20/10 

10/10 

12 June 

0/0 

40/20 

40/10 

0/0 

18 June 

0/0 

20/20 

20/20 

30/10 

24 June 

40/10 

10/0 

40/20 

10/0 

2 July 

— 

— 

20/0 

0/0 

15 July 

90/0 

60/10 

50/0 

0/0 

29 July 

100/60 

10/0 

60/20 

30/0 

12 August 

— 

— 

70/20 

10/0 

26 August 

90/70 

0/0 

40/0 

20/0 

19 September 

— 

— 

40/30 

10/0 

5 October 

0/0 

10/0 

60/20 

40/30 


LIN = Long Island Native. 

LI1 = Long Island Inbred. 

DBN = Delaware Bay Native. 

DBI = Delaware Bay Inbred. 

Stauber (1950), the Long Island Sound population was 
distinct because it began spawning when temperature 
reached 16.4°C, while other populations spawned at 
higher temperatures. In this study, both Long Island 
groups began spawning when water temperature was be¬ 
tween 15 and 20°C, while spawning in the Delaware Bay 
groups began when water temperature was about 25-28°C 
(Fig. 1). Thus there is little difference in the temperature 
at which Long Island oysters spawn in their native setting 
or in Delaware Bay, even after six generations of inbreed¬ 
ing in Delaware Bay. 

Gonadal development is an energy demanding process, 
with necessary nutrients coming either from stored re¬ 
serves, recently ingested food, or both (Sastry, 1979; Bar¬ 
ber and Blake, 1991). Site-specific variation in gonadal 
development of Afvtilus edulis, A. irradians , and Placo- 
pecleti magellanicus has been reported to be due to adap¬ 
tion (non-genetic) to local variations in environmental 
factors, most notably food availability (Newell et al. , 1982; 
Bricelj et al., 1987; MacDonald and Thompson, 1988). 
The fact that at the same site in Delaware Bay, both inbred 
and native Long Island oysters have distinctly different 
gonadal cycles than Delaware Bay oysters, effectively 
eliminates food abundance as the sole regulator of gonadal 
development in C. virginica . Perhaps gonadal develop¬ 
ment in oysters is initiated in the presence of adequate 
food supplies, at a genetically determined temperature. 


The results of this study and previous studies establish 
the existence of populations of C. virginica along the east 
coast of the LInited States that are genetically distinct based 
on differences in the timing of cycles of gonadal devel¬ 
opment and spawning. These intraspecific differences oc¬ 
cur between populations located north and south of Long 
Island, New York. Loosanoff and Nomejko (1951) found 
that oysters from Massachusetts and Long Island Sound 
spawned successfully in Milford, Connecticut, while oys¬ 
ters from New Jersey and Virginia did not spawn or only 
partially spawned. Similarly, Loosanoff (1969) saw a dif¬ 
ference in the temperature at which gonadal development 
was initiated between oysters from north of Long Island 
and south of Long Island. Oysters transplanted from 
Chesapeake Bay (Maryland) to Florida conformed to local 
spawning conditions after one year (Butler, 1955), indi¬ 
cating that observed differences in spawning temperatures 
between these locations were non-genetic. 

Genetic discontinuities along the distributional range 
of C. virginica have been reported at Brownsville, Texas, 
based on electrophoretic variation in proteins encoded by 
nuclear genes (Buroker, 1983) and along the mid-Atlantic 
coast of Florida based on restriction site variation in mi¬ 
tochondrial DNA (Reeb and Avise, 1990). It is noteworthy 
that populations of oysters exist that have genetically dis¬ 
tinct gonadal maturation cycles and yet are indistinguish¬ 
able based on both allozyme frequencies and mtDNA re¬ 
striction site variation. The fact that discrepancies exist 
between allozyme frequencies, mtDNA patterns, and 
physiological processes such as reproduction, point out 
that there is no simple means of delineating genetic 
boundaries among oyster populations (see review by 
Gaffney, 1991). 

The existence of stocks of oysters having genetically 
distinct reproductive cycles has implications for the oyster 
fishery. Oysters transplanted from a northern location 
(north of Long Island) to a southern location (south of 
Long Island), requiring lower temperatures to initiate go¬ 
nadal development, would spawn earlier than the local 
population. In the case of Long Island oysters in Delaware 
Bay, there is a potential one month overlap in spawning 
periods during which inter-breeding could occur. For the 
most part, however, the spawning periods would remain 
distinct. Oysters transplanted from southern to northern 
locations might not experience temperatures high enough 
to initiate gonadal development or spawning. In this case, 
oysters with unspawned or partially spawned and resorb¬ 
ing gametes would have little reproductive value and 
would be of poorer quality and lower commercial value. 
Hatchery production of oysters would be enhanced by the 
ability to obtain naturally conditioned broodstock over a 
wider time interval than what is available with local stocks 
alone. 
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A LONG. ISLAND NATIVE 


B LONG ISLAND INBRED 




C DELAWARE BA if NATIVE 


D DELAWARE BAY INBRED 



Figure 2. Mean gonadal area index (± 1 SD) for four groups of Crassostrea virginica: A = Long Island 
Native (LIN); B = Long Island Inbred (LI 1); C = Delaware Bay Native (DBN); D = Delaware Bay Inbred 
(DBI). Filled lines are spline interpolations. 


Table II 

Statistical comparison (l-test) of coefficients a and ft from the quadratic 
model for each possible combination of oyster, CrassosLrea virginica, 
groups H 0 : a x = a ) and ft t = ft } , where i and j are di fferent groups 
of oysters 


Groups 

1 J 

l 

(«, = Ofj) 

t 

(A = A) 

P* 

LIN—DBN 

6.17 

7.72 

P < 0.002 

LIN—LII 

1.39 

0.44 


LIN—DBI 

4.01 

6.33 

P < 0.02 

DBN—LII 

5.50 

6.66 

P < 0.005 

DBN—DBI 

2.62 

1.59 


LII—DBI 

2.89 

4.36 

P <s 0.05 


* If comparisons of both a and ft are rejected. 
LIN = Long Island Native. 

L1I = Long Island Inbred. 

DBN = Delaware Bay Native. 

DBI = Delaware Bay Inbred. 
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